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Abstract In this paper we study the multicast routing 
problem in all-optical WDM networks under the spare 
light splitting constraint. To implement a multicast ses- 
sion, several light-trees may have to be used due to the 
limited fanouts of network nodes. Although many mul- 
ticast routing algorithms have been proposed in order 
to reduce the total number of wavelength channels used 
(total cost) for a multicast session, the maximum num- 
ber of wavelengths required in one fiber link (link stress) 
and the end-to-end delay are two parameters which are 
not always taken into consideration. It is known that 
the shortest path tree results in the optimal end-to-end 
delay, but it can not be employed directly for multicast 
routing in sparse light splitting WDM networks. Hence, 
we propose a novel wavelength routing algorithm which 
tries to avoid the multicast incapable branching nodes 
(MIBs, branching nodes without splitting capability) in 
the shortest-path-based multicast tree to diminish the 
link stress. Good parts of the shortest-path-tree are re- 
tained by the algorithm to reduce the end-to-end delay. 
The algorithm consists of three steps: (1) a DijkstraPro 
algorithm with priority assignment and node adoption 
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is introduced to produce a shortest path tree with up 
to 38% fewer MIB nodes in the NSF topology and 46% 
fewer MIB nodes in the USA Longhaul topology, (2) 
critical articulation and deepest branch heuristics are 
used to process the MIB nodes, (3) a distance based 
light-tree reconnection algorithm is proposed to create 
the multicast light-trees. Extensive simulations demon- 
strate the algorithm's efficiency in terms of link stress 
and end-to-end delay. 

Keywords WDM Networks • Multicast Routing • 
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1 Introduction 

Multicast is a very efficient approach for one-to-many 
or many-to-many communication. A multicast session 
typically involves a source and a set of destinations. 
Typically in IP-based packet-switching networks a mul- 
ticast tree rooted at the source is constructed with 
branches spanning all destinations to accommodate a 
single source-based multicast session. In order to be 
able to multicast data in WDM optical networks, opti- 
cal switches need to have splitting capability. Note how- 
ever that optical switches with light splitters are always 
far more expensive to build than those without. Conse- 
quently when only a few nodes can support splitting the 
network is described as sparse light splitting [T]. Hence, 
multicast routing in WDM optical networks is very dif- 
ferent from that in IP-based packet switching networks 
and one must consider the constraint on node splitting 
capability in a practical optical network. To implement 
multicast in all-optical WDM networks, the light-tree 
concept was proposed in [2J. A light-tree is a point- 
to-multipoint generalization of a lightpath containing 
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one continuous lightpath from the source to each des- 
tination. To create light-trees, optical constraints must 
be respected. Without wavelength conversion the same 
wavelength must be used on all links in a light-tree. 
This is called the wavelength continuity constraint [3]- 
Moreover, two light-trees or lightpaths sharing a com- 
mon link must be assigned different wavelengths. This is 
known as the distinct- wavelength constraint [5] . Due to 
these physical constraints, supporting multicast routing 
in all-optical networks is a challenging task. 

For multicast routing in WDM optical networks many 
multicast light-tree computation algorithms have been 
proposed to reduce the total number of wavelength chan- 
nels used (i.e., the total cost), but the maximum num- 
ber of wavelengths required in one fiber link (i.e., the 
link stress) and the end-to-end delay are also important 
factors which should be taken into account. This is es- 
pecially true for time sensitive and bandwidth intensive 
multicast applications such as HDTV, VoIP and Video 
Conference. It is known that if a message is transmitted 
via the shortest path from the source to the destination, 
then, in general, the delay is minimal. Unfortunately, 
splitting nodes are very rare in optical networks due 
to their high cost and complex architecture. If a short- 
est path tree is directly applied for multicast routing in 
WDM networks, then there is a high probability that 
the branching nodes of this tree do not coincide with 
multicast capable switches. When this is the case dif- 
ferent wavelengths must be used to send messages from 
the source to different branches of a multicast incapable 
branching node, and the stress on the commonly used 
links will be very high. If the shortest paths are not 
used for all destinations, then a destination could find 
a longer path to the source (e.g., by connecting to a 
nearby splitting node), which implies a bigger end-to- 
end delay. Thus, a tradeoff must be found between link 
stress and end-to-end delay in order to obtain the best 
general performance. 

In this paper, a multicast routing algorithm con- 
sidering sparse light splitting is proposed which tries 
to avoid multicast incapable branching nodes in multi- 
cast light-trees. It aims to reduce both the link stress 
and the end-to-end delay. The significant aspects of 
this proposition are: (i) a DijkstraPro algorithm with 
priority assignment and node adoption, introduced to 
construct a shortest path tree with fewer multicast- 
incapable branching nodes, (ii) critical-articulation and 
deepest-branch heuristics are used to process the MIB 
nodes with the aim of reducing both link stress and end- 
to-end delay, (iii) a distance-based light-tree reconnec- 
tion algorithm is proposed to create a set of multicast 
light-trees with smaller end-to-end delay while keeping 
the same link stress and total cost. 



The rest of the paper is organized as follows: Re- 
lated work is reviewed in Section [5J The wavelength 
routing problem under the sparse light splitting model 
is formulated and some necessary definitions are given 
in Section |3l The multicast routing algorithm based 
on avoidance of multicast incapable branching nodes is 
proposed and simulated in Sections 2] and [5j Finally, 
a summary of results is made in Section [5] 

2 Related Work 

The difficulty of multicast routing in WDM networks 
with sparse light splitting has been addressed in many 
papers 4, 5, 6 , 7ll9l H0lH4lfT6] and various algorithms have 
been proposed. There are broadly three main categories 
according to the routing approaches they employ: Short- 
est Path Tree Based Routing (e.g., Reroute-to-Source 
and Reroute-to- Any [3]), Steiner-Based Routing (e.g., 
Member-Only [4: and Virtual-Source Capacity-Priority 
[3]) and Core-Based Routing (e.g., Virtual Source-based 
EE])- Essentially, the Shortest Path Based Routing ap- 
proach constructs the multicast tree by connecting the 
source to each destination individually using the ap- 
propriate shortest path in order to minimize the per- 
source-receiver path cost. The objective of the Steiner- 
Based Routing schemes, however, is to minimize the 
overall cost of the multicast light-trees. The Core-Based 
Routing algorithm connects a subset of nodes, called 
core nodes, which have both light-splitting and wavelength- 
conversion capacities. The multicast session is then es- 
tablished with the help of this core structure [3[7]. 
To the best of our knowledge from the literature [4], 
in WDM networks with sparse-splitting and without 
wavelength conversion the Member-Only algorithm yields 
the approximate minimal cost and the best link stress, 
while the Reroute-to-Source algorithm yields the opti- 
mal delay. 

In Reroute-to-Source, a multicast tree is first gener- 
ated to span all destinations, for example by comput- 
ing the shortest path tree with the Dijkstra algorithm. 
Then, it checks the light splitting capability for each 
branching node in the shortest path tree. If a branching 
node is a node with splitting capability, then no mod- 
ification is needed. But if it is a multicast incapable 
branching node (i.e., it has at least two direct children 
while it has no splitting capability), then only one di- 
rect child can be kept, which is chosen arbitrarily. All 
other direct children (sub-trees) must be connected to 
the source through the shortest path, each on a different 
wavelength. It is obvious that the end-to-end delay of 
Reroute-to-Source is minimal. The stress of the link is 
very high, however, because downstream branches of a 
multicast incapable branching node must be connected 
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to the source using the same shortest path on several 
different wavelengths. Note that there may actually be 
some longer paths leading to the source available on the 
same wavelength that are not used. 

In Reroute-to- Any, similarly to Reroute-to-source, a 
shortest path tree is first computed for all the destina- 
tions, and for each multicast-incapable branching node 
one downstream branch is kept and the others are cut. 
Finally, the cut destinations are reconnected to the mul- 
ticast light-tree via a multicast-capable node or a leaf 
multicast-incapable node in the light-tree if possible. If 
this is not possible they are reconnected to the source on 
different wavelengths. Although its link stress and total 
cost are better than the Reroute-to-Source and its av- 
erage end-to-end delay is superior to Member-Only, the 
algorithm is still not entirely satisfactory and should be 
improved to take traffic with QoS requirements into ac- 
count. It seems that no algorithm has been proposed to 
decide which branch of the multicast-incapable branch- 
ing nodes should be kept and what kind of reconnection 
algorithm can be used to reconnect the cut destinations. 

In Member-Only, each iteration adds the nearest 
destination to the multicast tree using the shortest path. 
The shortest path must not include any non-leaf mul- 
ticast incapable nodes in the light-tree under construc- 
tion. It is a modification of the well known Takahashi- 
Matsuyama heuristic of the Steiner problem [HI]- Al- 
though its total cost approaches the optimum, it is very 
likely that most of the destinations are connected to the 
source via a node far away from the source. As a result, 
its average end-to-end delay is high and the diameter 
of the multicast tree is always very large. 

3 Multicast Routing Under Sparse Light 
Splitting Constraint 

3.1 System Model and Problem Formulation 

One-to-many communication in all-optical WDM net- 
works is now examined. Network nodes equipped with 
light splitters are assumed to be sparse because of their 
complex architecture and expensive cost, with a pres- 
ence normally below 50% [13]. Furthermore, it is as- 
sumed that costly wavelength converters are not avail- 
able. With no loss of generality, the splitting capability 
of a multicast capable node is assumed to be infinite 
by supposing correct use of optical amplifiers [TS]. A 
spare light splitting WDM network can be modeled by 
an undirected graph G(V,E,c,d). Each node v € V is 
either a multicast incapable node (without splitting ca- 
pability) or a multicast capable node (equipped with 
light splitters) . Each edge e £ E is associated with two 
weight functions c(e), d(e). c(e) represents the cost of 



fiber link e, and d(e) denotes the propagation delay in 
fiber link e. Both of them are additive along a light- 
path LP(u,v). We consider the arrival of a multicast 
session ms(s,D), which requires a simultaneous com- 
munication from the source s to be set up to a group of 
destinations D. Due to the sparse-splitting constraint 
together with the wavelength-continuity constraint, one 
light-tree may not be sufficient to cover all destinations. 
Assume k light-trees LTi(s,Di) should be built for a 
multicast session ms{s,D) where i £ [l,k] and Di de- 
notes the set of destinations exclusively served in the 
i th light-tree. Since these k light-trees are not edge dis- 
joint, different wavelengths must be assigned for each 
light-tree. Thus, the number of wavelengths required for 
ms(s,D) (i.e., link stress) equals the number of light- 
trees built. 

Stress[ms(s, D)] — k (1) 

The total number of wavelength channels used (i.e., 
total cost) for ms(s, D) can be calculated as 

c[ms(s,D)]= E c ( e ) ( 2 ) 

i£[l,fc] e£LTi(s,Di) 

Nowadays multimedia services such as HDTV, VoIP, 
Video Conference and Video on Demand are widespread 
in the Internet. They are delay sensitive and bandwidth 
intensive. Consequently, the link stress and the delay 
are two important parameters for multicast light-tree 
selection in WDM optical networks. If the link stress is 
very high for one multicast session, fewer wavelengths 
can be allocated for other sessions. This may lead to a 
reduction of network throughput. 

Besides this, power loss and noise are two other 
challenging problems in all-optical networks. Although 
power loss can be compensated by appropriate place- 
ment of optical amplifiers in fibers and cross-connects, 
noise resulting from amplification can cascade and is 
hard to remove without electronic processing. It is prac- 
tical to limit the length of a path in order to decrease 
the number of amplifiers jlOj . In addition, optical net- 
works are becoming increasingly widespread in the In- 
ternet Backbone. Although optical messages are trans- 
mitted from the source to the destination at a very 
high speed, the nodes in WDM optical networks can 
be distributed around the world. When this is the case, 
the end-to-end delay will not be negligible and delay- 
sensitive traffic will require special treatment. Let LP(s, di) 
be the lightpath between the source s and the destina- 
tion di in the light-trees built for a multicast session 
m(s,D), the average end-to-end delay and the maxi- 
mum end-to-end delay can be defined as follow: 

AverDelay[ms(s,D)] = -Lx2 ]T d(e) (3) 
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MaxDelay[ms(s, D)] = max 



E 

e£LP(s,di 



d(e) 



(4) 



Note that end-to-end delay and link stress cannot 
be minimized simultaneously. If the shortest path tree 
is directly used for multicast routing, although its delay 
is optimal, its link stress is generally very high. When 
an approximated Sterner tree is employed to build the 
multicast light-trees (using the Member-Only algorithm 
[4]), the link stress is good, but the end-to-end delay is 
intolerable: an approach that produces a tradeoff solu- 
tion needs to be found. In order to minimize the end- 
to-end delay, the shortest path tree can be considered 
a good starting point for the construction of multicast 
light-trees. In order to improve the link stress, the num- 
ber of MIB nodes in the shortest path tree can be re- 
duced by making some destinations communicate with 
the source using longer paths. Putting this approach 
into practice, a multicast routing algorithm based on 
avoidance of MIB nodes is proposed for WDM net- 
works with sparse light splitting. To simplify the ob- 
jective metrics, the same edge cost and delay functions 
are applied throughout the network, and without loss 
of generality they are given as: 



c(e) = 1 unit cost d(e) = 1 unit delay 



3.2 Useful Definitions 



(5) 




Before describing our proposed multicast routing algo- 
rithm, some necessary definitions are introduced below. 

Definition 1: MI and MC nodes 

MI nodes: multicast incapable nodes are nodes which 
cannot split light signals, but have a Tap- and- Continue 
(TaC [H]) capability. That is to say, they can tap 
a small amount of optical power from a wavelength 
channel while forwarding it to only one output link. 
MC nodes: multicast capable nodes are nodes which 
are equipped with light splitters which can distribute 
the incoming message to all of the outgoing ports. 

In the illustration figures (except those for network 
topology, Figures [2] and [TU]) MI nodes are denoted 
by a rectangle while MC nodes are denoted by cir- 
cles. 

Definition 2: Multicast Incapable Branching Node (MIB 
node) 

MIB nodes have no splitting capability, but lead to 
several downstream branches in a multicast tree. An 
MIB node's out degree in the multicast tree is not 
less than two. Once an MIB node has forwarded the 
message to one branch it is incapable of forwarding 



Fig. 1 Critical Articulation Node 

it to another branch using the same wavelength. 

Definition 3: Set MC.SET, MLSET and D 

A set of light-trees may be required by a multicast 
session. For a multicast light-tree T under construc- 
tion, 

MC_SET: includes the MC nodes and the leaf MI 
nodes in T. They may be used to span T . Hence, 
nodes in MC_SET are also called connector nodes 
in T . 

MI_SET: includes only non-leaf MI nodes in T which 
are unable to connect a new destination to T. 
D: includes unvisited multicast members which are 
neither joined to the current multicast light-tree T 
nor to the previously constructed multicast light- 
trees. 

Definition 4: Constrained Path (CP) and Shortest Con- 
strained Path (SCP) 

A Constrained Path CP(u, T) between a node u and 
a tree T is defined as the shortest path SP(u,v) 
from node u to a connector node v in T, such that 
SP(u, v) does not traverse any node belonging to 
the MLSET of T. 

CP(u,T) = {SP(u,v)\v e MC.SET.and 

Vrr G SP(u, v),x& MLSET} 

Accordingly, the shortest of all possible Constrained 
Paths CPSet(u, T) is called the Shortest Constraint 
Path SCP{u,T). 



(6) 



c[SCP(u,T)} = min c[CP(u,T)] (7) 

CP(u,T)eCPSet(u,T) 

There may be several SCP(u, T) from u to T with 
different connector nodes v. 



Definition 5: Connection Constraint Node (CC node) 
and Critical Articulation Node (CAN) 
If node u is a CC node, there must be an interme- 
diate node which is included in all the paths from 
u to the source s. This intermediate node is called 
the critical articulation node: CAN(u, s). In other 
words, a CC node u cannot reach the source s with- 
out node CAN(u, s). 

For example, in Figure [TJ node CAN separates the 
network into two parts. Node d and source s are in 
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Algorithm 1 Avoidance of MIB Nodes for Multicast 
Routing 

Require: A multicast session ms(s,D) 
Ensure: A set of multicast light-trees for ms(s, D) 
1: Use the DijkstraPro algorithm to construct the shortest path 
tree SPT rooted at the source s. Prune all the non-destination 
leaf nodes and the nodes which do not lead to any destina- 
tions. 

2: Use the Deepest Branch and the Critical Articulation Heuris- 
tics to process the MIB nodes in SPT. 

3: Use the distance based light-trees reconnection algorithm to 
create the required set of light-trees for ms(s, D). 



different parts. Without node CAN, d is not able to 
communicate with s. So d is a CC node, and node 
CAN is the CAN(d, s). 

4 Avoidance of MIB Nodes for Multicast 
Routing 

The avoidance of MIB nodes multicast routing algo- 
rithm can be viewed as a post-processing of the short- 
est path tree (SPT). Due to the presence of MIB nodes 
in a shortest path tree a single wavelength may not 
be sufficient to cover all destinations and thus several 
wavelengths may be required to accommodate the mul- 
ticast group. Thus, MIB nodes should be avoided in 
order to decrease the link stress. If there are no MIB 
nodes in the shortest path tree, then the tree is an opti- 
mal multicast light-tree with both minimum end-to-end 
delay and minimum link stress. If this is not the case 
some processing must be done on the MIB nodes. The 
proposed algorithm consists of three main steps: the 
shortest path tree construction step, the MIB nodes 
processing step and the multicast light-tree reconstruc- 
tion step. In the first step, an enhanced DijkstraPro 
algorithm making use of the priority method and node 
adoption is introduced to construct a multicast short- 
est path tree with fewer MIB nodes and smaller link 
stress. In the second step the MIB nodes in the short- 
est path tree are processed: deepest branch and critical 
articulation heuristics are proposed to keep only one 
downstream branch of MIB nodes in an attempt to re- 
duce both the link stress and the end-to-end delay. In 
the last step the distance-based light-tree reconnection 
algorithm (which can also reduce end-to-end delay) is 
applied to create the multicast light-trees. 

4.1 Construction of SPT and DijkstraPro Algorithm 

First of all, a shortest path tree rooted at the source is 
constructed for all nodes in the network. Then, accord- 
ing to the multicast session, non-destination nodes and 



nodes that do not lead to any destination are pruned 
from this shortest path tree. 

Generally, Dijkstra's algorithm is employed to build 
the shortest path tree. In the Dijkstra algorithm, a node 
is said to be labeled permanently [11] if its shortest path 
to the source has been found. Otherwise it is said to be 
tentatively labeled [TT]. Initially, only the source s is 
permanently labeled and all the other nodes are ten- 
tatively labeled. In each iteration, the node with the 
shortest distance to the source among all the tenta- 
tively labeled nodes is chosen and labeled permanently. 
It is worth noting that in one iteration there may be 
several nodes that have the same shortest distance to 
the source, here we call them candidate nodes and the 
distance is referred to as their level. However, according 
to the Dijkstra algorithm we should label only one of 
the candidate nodes permanently in order to update the 
distances of the other nodes. The question then, is how 
to choose the permanently labelled candidate node? In 
the traditional Dijkstra algorithm, it is chosen arbitrar- 
ily. But consider this situation: there are two candidate 
nodes at the same level; one is a MI node and another 
is a MC node; they share the same two adjacent nodes. 
If the MI candidate node is the first to be selected for 
permanent labelling then the two adjacent nodes will 
update their distances to the source, and thus will be 
connected to the source via this MI candidate node. The 
problem is that the MI candidate node cannot split the 
incoming signal to more than one outgoing port. As a 
result, it will become a MIB node in the shortest path 
tree. Alternatively, if the MC candidate node is the first 
to be permanently labeled then when the two adjacent 
nodes update their distances to the source they will 
be connected to it via this MC candidate node. Subse- 
quently, the MI candidate node is chosen to be perma- 
nently labeled. At this point, no adjacent node needs 
to update its distance and no adjacent node is left to 
be connected to the source via this MI candidate node. 
So, the risk that an MI candidate node will become an 
MIB node is reduced or even avoided. 

Due to the constraint on splitting capability, the 
traditional Dijkstra algorithm may not yield a favor- 
able result, but it can be improved with some modi- 
fication. Hence, the DijkstraPro algorithm with prior- 
ity and node adoption is presented. When building a 
shortest path tree using Dijkstra with several candi- 
date nodes at the same level the following operations 
are proposed: 

— Giving Higher Priority to MC Candidate Node 
The candidate node with multicast splitting capa- 
bility (MC candidate node) should be given higher 
priority than the MI candidate nodes due to the 
fact that they can connect many destination nodes 



to the tree without producing a MIB node. In other 
words, the probability that an MI candidate node 
will be used to connect more than one destination 
to the tree in latter iterations is greatly reduced. 
Refer to the NSF network in Figure [2] Nodes 1, 8 
and 10 are assumed to be MC nodes. A multicast 
session arrives: mi — {source: 10 | members: 1 ~ 
14}. If the Dijkstra algorithm is used then we can 
get the shortest path tree in Figure [3] There are 2 
MIB nodes in this shortest path tree. We can see, 
however, that nodes 1, 6, 7, 9 and 13 have the same 
shortest distance to the source node 10. So, they 
can be viewed as candidate nodes at the same level. 
And, if node 1 (an MC node) is promoted to a higher 
priority and chosen first to be permanently labeled, 
followed by 7, 9, 13 and 6, then we can get the new 
shortest path tree of Figure 2] which has only one 
MIB node. 

Giving High Priority to MI Candidate Nodes with 
Smaller Degree 

If there are no MC candidate nodes, then the candi- 
date node with the smaller degree is given a higher 
priority. This is because the probability that an MI 
candidate node with a smaller degree is a MIB node 
is very low (especially for candidate nodes with a 
degree of two) . That is to say, the number of nodes 
that remain to be connected to the source through 
other MI candidate nodes with higher degrees is very 
small. Consequently, the probability for a candidate 
node with a higher degree to become an MIB node 
is reduced. So, the average probability for a node to 
be a MIB node is slightly reduced. 




Fig. 2 NSF Network Topology 
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Fig. 3 The SPT for m\ constructed by the Dijkstra algorithm 




Fig. 4 The SPT for mi constructed by offering higher priority 
to MC candidate nodes 



Node Adoption 

At the stage when all candidate nodes at the same 
level have been permanently labeled, the following 
situation may occur: some MI candidate nodes con- 
nect only two direct children to the tree (i.e., MIB 
candidate nodes) while some candidate nodes are 
leaf nodes in the created tree. Thus, the possibility 
arises for a leaf candidate node to adopt one child 
from an MIB candidate node at the same level when 
the child can reach the source through the leaf can- 
didate node also. By doing this the creation of an 
MIB node can be avoided. Node adoption between 
the candidate nodes at the same level can result 
in a greatly reduced number of MIB nodes in the 
shortest path tree or in the balancing of the load of 
an MIB node. Typically a destination node should 
be given a higher priority when determining which 
nodes may be adopted. 



Fig. 5 The SPT after Node Adoption from Figure [4] 

Refer to the example in Figure @] It is obvious that 
nodes 11, 12 and 14 have the same least distance 
to the source node 10, hence they can be viewed as 
candidate nodes at the same level. After all of them 
have been permanently labeled we can see that node 
12 is an MIB node and node 14 is a leaf node. Note 
that nodes 13 or 9 can reach the source node 10 by 
the shortest path through both of nodes 12 and 14. 
Thus, one of them could be adopted by node 14, 
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and a new shortest path tree without an MIB node 
can be obtained as in Figure El 

4.2 Processing of the MIB nodes 

Due to the fact that the MIB nodes in the shortest path 
tree can forward the incoming message to one and only 
one outgoing branch, the existence of MIB nodes is the 
most important cause of high link stress. Thus, they 
should be processed and avoided. In the Reroute-to- 
Source algorithm [4] , all downstream links of MIB nodes 
are connected to the source through the reverse short- 
est path on different wavelengths which results in high 
link stress. Although the Reroute-to- Any algorithm is 
also proposed in the literature 0], there is no descrip- 
tion of how to keep one branch when processing the 
MIB nodes. So, in this paper, the deepest branch and 
critical articulation heuristics are employed to decide 
which branch should be kept in order to decrease the 
link stress and the end-to-end delay. 

4.2.1 MIB Pro 

— Critical Articulation Heuristic 

A CC node u can only communicate with the source 
through its CAN(u, s). In a multicast tree, if the 
CAN(u, s) is (unfortunately) an MIB node, then 
the branch containing u should be assigned a higher 
priority and kept when processing this MIB node. 
This is because there is no alternative path for u to 
reach the source without traversing its CAN(u, s). 
However, destinations in the other branches may 
find another path to the source which will not tra- 
verse this MIB node. In fact, CC and CAN(CC, s) 
nodes are very rare in real optical networks. How- 
ever, in the case that some nodes in the network 
have failed they may exist, and this heuristic will 
be very practical. In the network of Figure H2 node 
g?i is a CC node. The shortest path tree for multi- 
cast session rri2 = {source: s | destinations: d\ ~ 
} is given by Figure |U We can see that CAN(di, s) 
is an MIB node in the shortest path tree built for 
rri2 as plotted in Figure El hence it should be pro- 
cessed. If node di is disconnected from CAN(d\, s) 
and the branch leading to node di and d^ is kept, 
then two light-trees on two different wavelengths wq 
and wi are required as shown in Figure [3 But if the 
CC node d\ is kept and the other one is cut, then 
only one light-tree (or one wavelength) is needed as 
shown in Figure |8l 

— Deepest Branch Heuristic 

The deepest branch should also be assigned a higher 



Algorithm 2 Processing of MIB nodes Using Critical 
Articulation and Deepest Branch Heuristics 

1: Search all the MIB nodes in the shortest path tree 
2: for each MIB node do 

3: if No downstream branch contains a CC node then 
4: Keep the deepest branch 

5: else if Only one downstream branch contains a CC node 

& MIB node = CAN{CC, s) then 
6: Keep the branch with the CC node 

7: else if Several downstream branches contain CC nodes 

& MIB node = CAN(Cd , s), i = 1, 2, . . . then 
8: Keep the deepest branch with a CC node 

9: end if 
10: end for 

11: Delete the downstream branches of MIB nodes which are not 
kept 



□ --MI Q-MC 




Fig. 6 An example network with CAN nodes 



>.»„ It', 

LT, LT, 




Fig. 7 A Shortest Path Tree for 7712 




Fig. 8 Processing MIB nodes using the Critical Articulation 
Heuristic 

priority. This is because a destination far away from 
the source has difficulty finding a path to the source 
without traversing a non-leaf MI node in the tree. 
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Furthermore it is desirable to minimise the average 
end-to-end delay for a destination node far away 
from the source by choosing the shortest path to 
the source. To implement this step a breadth-first 
traversal algorithm can be employed. The worst case 
time complexity of this heuristic is O(N), where 
N = \V\ is the number of nodes in the network. 

4.2.2 MIBPro2 

In addition another method is proposed to process MIB 
nodes in the shortest path tree. This method deletes 
all the downstream branches of an MIB nodes without 
employing any heuristic. These two methods will be 
compared in Section [5] 



4.3 Reconnection of Multicast Light-trees 

After the MIB node processing step, the shortest path 
tree is divided into a disconnected forest containing a 
subtree plus several separated destinations. This dis- 
connected forest must be reconnected in order to ac- 
commodate all the multicast members. A Member-Only- 
like [1] light-tree connection method would be a good 
candidate to reconnect the multicast forest. The Member- 
Only algorithm always adds the destination nearest to 
the multicast light-tree using the shortest path, but this 
shortest path will not use any non-leaf MI node in the 
light-tree. In other words, at each iteration only the 
destination with the shortest SCP is connected to the 
light-tree using this SCP. As demonstrated in [4], the 
Member-Only algorithm can achieve the best link stress 
and the minimum cost, although its end-to-end delay is 
very large. It is worth noting that some improvements 
can be made to this algorithm to reduce the end-to-end 
delay to some extent while obtaining the same cost and 
the same link stress. The example below demonstrates 
how end-to-end delay can be reduced. 

A multicast session 777,3 = {source: 10 | destinations: 
6, 11, 13, 14} is required in the NSF network, Fig- 
ure [2] We assume that the first tree only contains the 
source node 10. According to the previously described 
member-only-like light-tree reconnection approach, the 
destination with the shortest SCP should be added to 
this tree first. The shortest paths for node 11 and node 
14 to the source have length 1. Without loss of general- 
ity, let us suppose node 14 is the first to be connected. 
Then, on the new tree, we can see that both SCPs for 
nodes 11 and 13 have the same length. Also without 
loss of generality, suppose node 13 is then connected. 
After that node 6 is chosen, and finally node 11. Fol- 
lowing these steps, the resultant multicast tree is given 
in Figure [HJa) . 



Algorithm 3 Distance Based Light-tree Reconnec- 
tion Algorithm 



10: 
11: 

12: 

13 
14 
15 

16: 

17 
18 
19 
20 
21 
22 
23 
21 
25 
26 
27 



T <— subtree obtained after MIB process 
MCSET <r- {MC nodes and leaf MI nodes in T} 
MI. SET <- {non-leaf MI nodes in T} 
D <— {destinations not in T} 
while (D j= <P ) do 
repeat 

Find the closest destination d £ D to T, such that 
its shortest path to T does not traverse any node in 
MLSET 

if there are several destinations satisfying equa- 
tion [8] then 

Select the destination nearest to s in network 

G as d 
end if 

if there are several connector nodes for d in 

MCSET satisfing equation [9] then 

Select the connector node nearest to s in T as 
c and choose the corresponding SCP 

end if 

T ^TUSCP(d,c) 

MCSET <r- MCSET U {d and MC- 
nodes on SP(d,c)} 

MLSET +- MLSET U {non-leaf MI- 
nodes on SP(d, c)} 
D <- D\d 

if c is an MI node then 

MCSET <- MCSET \ c 

MLSET «- MLSET U {c} 
end if 

until no destination can be added to T 
return T 

Begin a new tree T <— {s} 
MCSET <- {s} 
MLSET <- <j> 
end while 



dist{SCP(d, T)} -- 
dist{SCP(d,T)} ; 



min dist[SCP(di,T)] 
- dist{SP(d, connectovi)}, 



1,2, 



(8) 
(9) 





(b) 



Fig. 9 Two strategics for the reconnection of light-trees 

It is immediately apparent that node 11 can be con- 
nected to the tree via node 10 or node 6. Why do we 
not connect it through node 10 as in Figure |5fb)? The 
difference is that the connector nodes have different dis- 
tances to the source in the tree (for node 10 the distance 
is while that for node 6 is 3). In addition, it is even 
more interesting to consider Figure M, c ) ■ All three of 
these multicast trees have the same cost of 4 while hav- 
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Fig. 10 USA Longhaul Network 

ing different average delays: 10/4, 7/4 and 6/4. It is 
also simple to determine that following the addition of 
node 14 to the tree, if node 11 is added before node 13 
we can get the result in Figure [^c). 

So, from this simple example we have two strategies 
that reduce the average delay while maintaining the 
same cost and the same link stress. The distance based 
reconnection algorithm is developed from these obser- 
vations. If there are several nodes whose SCPs to the 
multicast tree have the same length, then these nodes 
should be added in the order of their distance to the 
source (the distance in the network): the nearer, the 
earlier. When the destination with the shortest SCP 
has at least two connector nodes in the subtree, it is 
better to use the connector node nearest to the source 
(the distance in the multicast light-tree under construc- 
tion), otherwise its end-to-end delay will be too large. 

5 Performance Evaluation and Simulation 

To ensure the effectiveness of our proposed multicast 
routing algorithm, two different network topologies are 
employed as test beds for the simulation: the 14 node 
NSF network in Figure [2] and the 28 node USA Long- 
haul network in Figure [TO] The fact that these net- 
works have been used as reference topologies in many 
papers [51l71 [TUl[T4lHo1[T7] is the reason for their selec- 
tion. 

5.1 Performance of the DijkstraPro Algorithm 

To demonstrate the superiority of the DijkstraPro al- 
gorithm it is compared with the traditional Dijkstra 
algorithm using the following two parameters: 

— N: the number of the MIB nodes in the shortest 
path tree. 



Table 1 Comparision of Dijkstra and DijkstraPro algorithms in 
the NSF Network in Figure [2] 
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source 




MC 
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DijkstraPro 
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1 


2 


4 
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3 
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2 


4 


3 


3 


2 
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4 
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3 


3 


2 


1 


2 
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3 


2 


3 


2 


3 


2 


3 


2 


7 


5 


5 


4 


4 


3 


3 


2 


2 


8 


4 


4 


2 


3 


3 


2 


1 


2 


9 


4 


3 


3 


3 


3 


3 


2 


3 


10 


3 


2 


1 


2 


2 


2 


1 


2 


11 


3 


4 


3 


4 


1 


2 


1 


2 


12 


4 


3 


4 


3 


2 


2 


1 


2 


13 


3 


4 


2 


4 


2 


2 


1 


2 


14 


4 


3 


4 


3 


2 


2 


1 


2 


Average 


3.64 


3.43 


2.79 3.07 


2.43 2.29 


1.5 2.14 



— S: the maximum number of wavelengths required in 
one fiber link to cover all destinations in the shortest 
path tree (i.e., link stress of the SPT). 

In each comparison, two conditions are considered. 
Condition 1 only regards the source to be an MC node, 
while Condition 2 regards nodes with a high degree to 
be MC nodes. The reason for choosing these two con- 
ditions can be explained as follows. In Condition 1, as 
only the source is an MC node, MC candidate node 
priority is not applied. Thus, the result in Condition 
1 checks the merit of the node adoption operation in 
the DijkstraPro algorithm. As stated in |18l[191 , one ap- 
proach could be to place splitters at the nodes with high 
degree. Thus, nodes with high degree are treated as MC 
nodes in Condition 2. In this condition, MC candidate 
node priority is applied, and the overall performance 
can be verified. 

In Table [TJ we evaluate the performance of 14 short- 
est path trees rooted at each node of the NSF network. 
Source ID denotes the root of the shortest path tree 
built. Two conditions are considered: 

Condition 1 (only the source is an MC node) 

The average number of MIB nodes in the shortest 
path tree constructed by the DijkstraPro algorithm 
is 0.85 less (23%) than when applying the tradi- 
tional Dijkstra algorithm and the link stress is 0.36 
smaller. This result confirms that the node adoption 
operation in the DijkstraPro algorithm is effective. 

Condition 2 (nodes 6, 10 and the source are MC nodes) 
In the NSF network, node 6 and node 10 have a 
high degree of 4), so they can be assumed to be MC 
nodes which are very useful for multicast sessions. 
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The DijkstraPro algorithm produces a shortest path 
tree with fewer MIB nodes and smaller link stress 
for this condition also. The average number of MIB 
nodes is 0.93 less (38%) and the link stress is 0.15 
smaller. 

In Table [2] we also provide the performance of 28 
shortest path trees rooted at each node in the USA 
Longhaul network. 

Condition 1 (only the source is an MC node) 

The DijkstraPro algorithm results in 1.75 (29%) 
fewer MIB nodes on average than the traditional Di- 
jkstra algorithm, and the link stress of the shortest 
path tree built by DijkstraPro is 1.64 smaller. This 
signifies that the effectiveness of the node adoption 
operation is independent of network topology. 

Condition 2 (nodes 10, 12-15, 18, 21, 26 and the 
source are MC nodes) 

In the USA Longhaul network, nodes 10, 12—15, 18, 
21 and 26 have a degree equal to or above 4, so they 
are regarded as the MC nodes in this condition, the 
DijkstraPro algorithm can also produce a shortest 
path tree with fewer MIB nodes and smaller link 
stress. The average number of MIB nodes is 0.78 less 
(46%) and the average link stress is 0.43 smaller. 

Moreover, it is evident that when all the nodes in a 
WDM network are MC nodes, none of the shortest path 
trees constructed by the Dijkstra or the DijkstraPro al- 
gorithm will have any MIB nodes and their link stress 
will always be 1. So, it is obvious that when the ratio of 
MC nodes in the network is very high the improvement 
to be gained by using the DijkstraPro algorithm is not 
significant. But when the MC nodes are very sparse its 
performance is much better than the traditional Dijk- 
stra algorithm, not only in terms of the number of MIB 
nodes but also in terms of the link stress. This justifies 
our introduction of the DijkstraPro algorithm in the 
shortest path tree construction step for the implemen- 
tation of our proposed multicast routing algorithm. 

5.2 Performance of the Avoidance of MIB Nodes 
Based Multicast Routing Algorithm 

There is no mention in the literature for the Reroute- 
to- Any |3] algorithm of a technique to determine which 
branch of a MIB node should be cut, and which algo- 
rithm should be used to reconnect the cut destinations. 
In our simulation an arbitrary branch is assumed to be 
kept and a Member-Only-like [4] reconnection method 
is employed. 

To evaluate the performance of the proposed mul- 
ticast routing algorithm based on avoidance of MIB 



Table 2 Comparision of Dijkstra and DijkstraPro algorithms in 
the USA Longhual Network of Figure [TO] 
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nodes (MIBPro/MIBPro2), the following four metrics 
are used to measure the quality of the multicast light- 
trees built for a multicast session: 

— link stress 

— average end-to-end delay 

— maximum end-to-end delay 

— total cost 

In addition, each multicast session has only a single 
source. Each network node is selected as the source of 
a multicast session in turn. The destinations of a multi- 
cast group are distributed independently and uniformly 
through the network. For a given source and a given 
multicast group size, 100 random multicast sessions are 
generated. Hence, the result of each point in the simu- 
lation figures is the average of 100 x |V| computations. 
In addition, Reroute-to-Source (R2S), Reroute-to- Any 
(R2A) and Member-Only (MO) are also implemented 
for comparison. 



11 




Fig. 11 Link Stress vs. Multicast Group Size in (a)NSF Network (b)USA Longhaul Network 




Fig. 12 Average End-to-End Delay vs. Multicast Group Size in (a)NSF Network (b)USA Longhaul Network 




Fig. 13 Total Cost vs. Multicast Group Size in (a)NSF Network (b)USA Longhaul Network 




Fig. 14 Maximum End-to-End Delay vs. Multicast Group Size in (a)NSF Network (b)USA Longhaul Network 
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5.2.1 Effect of Group Size (Number of Multicast 
Members ) 

Here we study the performance of the proposed algo- 
rithm versus multicast group size. As mentioned in sub- 
section 15. 1[ nodes with high degree have a high prob- 
ability of being MC nodes. To simplify the simulation 
in this part, we regard these nodes as MC nodes and 
only change the group size to evaluate the quality of 
light-trees built by MIBPro/MIBPro2 multicast rout- 
ing algorithms. 

In the NSF network, nodes 6, 10, and the source are 
set as MC nodes. The simulation results in the NSF 
network are plotted in Figures ITTIfCT b). As shown in 
Figure fTTT a). we can see that when the group size is 
above four, MIBPro achieves better link stress than 
R2A. The link stress of MIBPro2 is also much smaller 
than MIBPro. Figures fl2Ta) and HITa) show that the 
average end-to-end delay and maximum end-to-end de- 
lay of MIBPro is only second to the optimal result of 
R2S. As multicast group size grows the improvement 
of end-to-end delay returned by MIBPro compared to 
R2A becomes more and more significant. Moreover, while 
the total costs of R2A, MIBPro and MIBPro2 are al- 
most the same, R2S results in the highest and MO re- 
sults in the lowest total cost. 

In the USA Longhaul network, nodes 10, 12~15, 
18, 21, 26, and the source are set as MC nodes. Fig- 
ures njHUb) have compared the performance of those 
five algorithms in this topology. The link stress of the 
five algorithms are almost the same and very near to 
1 according to Figure QjTb). This is because the ratio 
of MC nodes is very high (32%) in this configuration. 
The end-to-end delay for the MIBPro algorithm is very 
close to the optimum (R2S). To our surprise, MIBPro 
obtains almost the same maximum end-to-end delay as 
R2S. From the point of view of the total cost, R2A, 
MIBPro and MIBPro2 return the same value, which is 
the same outcome as the NSF network example. 

In both topologies the performance of R2S in terms 
of link stress and total cost is always the worst, while its 
performance in end-to-end delay is the best. Conversely 
the MO algorithm can achieve very good link stress and 
total cost, while its end-to-end delay is too large. 

From the simulation results above it can be seen 
that the MIBPro algorithm can provide nearly the same 
or even slightly better link stress than R2A. Its reduc- 
tion in average and maximum end-to-end delay com- 
pared to R2A becomes more obvious when the group 
size is large. This is because the MC node priority 
mechanism, node adoption and distance based recon- 
nection do not affect the result when the group size 
is too small. Only when there are enough destinations 



can these strategies work well. Overall, however, the 
MIBPro algorithm achieves a good tradeoff between 
link stress and end-to-end delay. 

5.2.2 Effect of Splitting Capability (Number of MC 
nodes ) 

The performance when the number of MC nodes varies 
have also been studied. According to the results of the 
previous section, MIBPro is more advantageous when 
the multicast group size is large. Thus, the multicast 
group size is set at a large value while only the num- 
ber of MC nodes is changed in the simulation of this 
part. The MC nodes are assumed to be independently 
and uniformly distributed in the topology. The multi- 
cast group size is set to 12 in the 14 nodes NSF network 
and set to 21 in the 28 nodes USA Longhaul network. 
The numeric results are plotted in Figures [TSTfTSl Ac- 
cording to these figures, when MC nodes are sparse, 
(1) MIBPro achieves much better performance in terms 
of link stress, average end-to-end delay and maximum 
end-to-end delay relative to R2A while producing the 
same cost as R2S. (2) MIBPro2 results in both lower 
link stress and total cost than R2A. Its link stress is 
even better than MO in the Longhaul network. How- 
ever, its end-to-end delay is either better or worse than 
R2A. 

These results indicate that our proposed MIBPro 
algorithm works well in the case of sparse light splitting. 
When the ratio of MC nodes is large, there are fewer 
MIB nodes in the shortest path tree and as a result 
MIBPro's advantage is less significant. 

6 Conclusion 

Due to physical constraints multicast routing in WDM 
networks with sparse light splitting is not easy. Many 
multicast routing algorithms have been proposed to find 
the Steiner based light-tree in WDM networks with the 
minimum cost, but it has been proved that this prob- 
lem is NP-hard. Applications with QoS requirements 
are becoming more and more popular in the Internet. 
The bandwidth (or number of wavelengths supported 
per fiber link in WDM networks) and the end-to-end 
delay are two important parameters for QoS. Hence, a 
multicast routing algorithm based on avoidance of MIB 
nodes is presented for traffic with QoS requirements in 
WDM networks in order to decrease the link stress and 
the end-to-end delay. The algorithm retains the parts 
of the shortest path tree which result in optimal end- 
to-end delay for at least some multicast members. In 
order to reduce the number of MIB nodes and the link 
stress in the construction of the shortest path tree step 
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Fig. 15 Link Stress against the Number of MC Nodes in (a)NSF Network (b)USA Longhaul Network 




Fig. 16 Average End-to-End Delay against the Number of MC Nodes in (a)NSF Network (b)USA Longhaul Network 




Fig. 17 Total Cost against the Number of MC Nodes in (a)NSF Network (b)USA Longhaul Network 




Fig. 18 Maximum End-to-End Delay against the Number of MC Nodes in (a)NSF Network (b)USA Longhaul Network 
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the DijkstraPro algorithm is presented, where a higher 
priority is assigned to MC candidate nodes and node 
adoption is performed between the candidate nodes at 
the same level. To keep one branch of MIB nodes in 
the shortest path tree, critical articulation and deepest 
branch heuristics are introduced. Finally, the distance- 
based light-tree reconnection algorithm is developed to 
rejoin the multicast light-forest. The first part of the 
simulation in Section[5]shows that the DijkstraPro algo- 
rithm is a better tool for shortest path tree construction 
in all-optical networks than the traditional Dijkstra al- 
gorithm. It can really reduce the number of MIB nodes 
and the link stress of the shortest path tree. Moreover, 
the second part of the simulation proves that the pro- 
posed MIBPro algorithm yields good performance in 
terms of link stress when MC nodes are very sparse. 
In addition, when the group size is large enough it is 
able to improve the average and maximum end-to-end 
delay dramatically giving a result very close to the op- 
timal Reroute-to-Source algorithm solution [3] . To sum 
up, the proposed algorithm is a good tradeoff between 
link stress and end-to-end delay for multicast routing 
in sparse light splitting WDM networks. 
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